! ******** û ***************
/clear
/filn, E132 
/title, Fatigue analysis of upper head of adsorber
/units,si                            !ùʵλ
!   
Di=4400                              ! 豸ھ
tc=16                                ! Ͳ弰ͷں
Hc=5*nint(sqrt(Di/2*tc)/10)*10       ! ģͲ峤ȣͷֱ߸߶)

dn1=630                   ! ӹ1⾶
dn2=530                   ! ӹ2⾶
dn3=426                   ! ӹ3⾶

tn1=8                     ! ӹ1ں
tn2=8                     ! ӹ2ں
tn3=8                     ! ӹ3ں

ln1=350                   ! ӹ1߶  
ln2=250                   ! ӹ2߶ 
ln3=200                   ! ӹ3߶

x0=1000                   ! ӹ1,2λ

nt=2                      ! ȷʷ
nme=30                    ! Բͷʷ
nh=20                     ! ӹܻʷ

p1=-0.1                   ! ѹ
p2=0.25                   ! ѹ
pax11=-p1*dn1**2/((dn1+2*tn1)**2-dn1**2)         ! ѹ½ӹ1˲ƽغ
pax12=-p1*dn2**2/((dn2+2*tn2)**2-dn2**2)         ! ѹ½ӹ2˲ƽغ
pax13=-p1*dn3**2/((dn3+2*tn3)**2-dn3**2)         ! ѹ½ӹ3˲ƽغ
pax21=-p2*dn1**2/((dn1+2*tn1)**2-dn1**2)         ! ѹ½ӹ1˲ƽغ
pax22=-p2*dn2**2/((dn2+2*tn2)**2-dn2**2)         ! ѹ½ӹ2˲ƽغ 
pax23=-p2*dn3**2/((dn3+2*tn3)**2-dn3**2)         ! ѹ½ӹ3˲ƽغ
! ********* ǰ  ************
/PREP7 
et,1,95                   ! 嵥Ԫ
mp,ex,1,2e5               ! ϵĵģ
mp,nuxy,1,0.3             ! ϵĲɱ

! ********* ģ  ************
CSWPLA,11,1,0.5,1,         ! Բϵ
K,1,-Di/2,,,               ! Բһ˵
K,2,Di/2,,,                ! Բһ˵
L,1,2                      ! Բ
csys,0                     ! ֱϵ
k,3,kx(2)+tc,,             ! ͷںؼ
L,2,3                      ! ɷͷں
ADRAG,2, , , , , ,1        ! ɷͷںԲԲ
BLC4,dn3/2-tn3, ,tn3,Di/4+tc+ln3,   ! Ľӹܽ  
APTN,all                            ! 򻥷
adele,4,5,,1                        ! ɾ

*GET,KPMAX,KP,,NUM,MAX              ! ȡڵ
k,KPMAX+1,0,0,0                     ! һؼ
k,KPMAX+2,0,hc,0                    ! һؼ
VROTAT,ALL,,,,,,KPMAX+1,KPMAX+2     ! תԲͷĽӹ
NUMMRG,ALL, , , ,LOW                ! ϲ
NUMCMP,ALL                          ! ѹ

asel,s,loc,y,0                      ! ѡԲͷ
*GET,KPMAX,KP,,NUM,MAX              ! ȡڵ
k,KPMAX+1,0,0,0                     ! Ͳһؼ
k,KPMAX+2,0,-hc,0                   ! Ͳһؼ
l,KPMAX+1,KPMAX+2                   ! Ͳ  
VDRAG,all, , , , , ,74              ! Ͳ 
ldele,1,,,1                         ! ɾ

alls                                ! ȫѡ
wprot,,-90                          ! תƽ
wprot,45                            ! תƽ

*AFUN,deg                           ! 趨ǶȵλΪ
ang1=acos(x0/Di*2)                  ! ӹλý
fai=atan(1/2/tan(ang1))             ! ӹλý
y0=Di/4*sin(ang1)                   ! ӹλ
y0c=y0-x0/tan(fai)                  ! ӹ߽
l0c=x0/sin(fai)                     ! ӹ߽ľ
wpoff,,,y0c                         ! ƽƹƽ
wprot,,,-fai                        ! תƽ
CYL4, , ,dn1/2-tn1, ,dn1/2, ,l0c+tc+ln1     ! ɽӹ1

wprot,,,fai                                 ! תƽ
wprot,,-fai,                                ! תƽ
wprot,90                                    ! תƽ
CYL4, , ,dn2/2-tn2, ,dn2/2, ,l0c+tc+ln2     ! ɽӹ2
wprot,-90                                   ! תƽ

wprot,,fai,                                 ! תƽ
wpoff,,,-y0c+Di/10                          ! ƽƹƽ
vsel,s,,,17,18,1                            ! ѡӹ12
vsbw,all                                    ! ùиӹ12
vsel,r,loc,y,-hc*10,Di/10                   ! ѡӹ12ಿ
vdele,all,,,1                               ! ɾӹ12ಿ

vsel,s,loc,y,0,hc*10                        ! ѡӹ12ͷز                
vsel,r,loc,x,-Di,0                           
VPTN,all                                    ! ӹ12ͷ廥
vdele,23,24,,1                              ! ɾӹ12ಿ
vdele,17,18,,1                              ! ɾͷಿ  

alls                                        ! ȫѡ
vsel,s,loc,x,-Di,0                          ! ѡӹ1Ľӹܡͷز  
vsel,r,loc,z,0,Di                           
wprot,,90,                                  ! תƽ
vsbw,all                                    ! иӹ1Ľӹܡͷز

vsel,s,loc,x,-Di,0                          ! ѡӹ2Ľӹܡͷز
vsel,r,loc,z,-Di,0                         
wprot,,,90,                                 ! תƽ
vsbw,all                                    ! иӹ2Ľӹܡͷز

lsel,s,length,,0,tc                         ! ѡȷ߶
LESIZE,all, , ,nt,, , , ,1                  ! 趨ȷʷ  
vsel,s,loc,y,-hc,0                          ! ѡͲ
vsweep,all                                  ! ʷͲ

vsel,s,loc,x,0,Di                           ! ѡxĲַͷĽӹ
vsel,u,loc,y,-hc,0                          ! ȥͲ                 
local,11,1,,,,,90                           ! ֲ
aslv                                        ! ѡ
lsla                                        ! ѡ
lsel,r,loc,x,dn3/2,Di                       ! ѡԲͷ
lsel,u,loc,x,dn3/2                          ! ȥӹ߶
lsel,u,loc,z,0                              ! ȥͲ߶
LESIZE,all, , ,nme, , , , ,1                ! 趨Բͷʷ                 
vsweep,all                                  ! ɨʷxĲַͷĽӹ    

alls                                        ! ȫѡ
lsel,s,loc,y,180                            ! ѡ180λ
lsel,r,loc,x,dn3/2,Di                       ! ѡĽӹ
lsel,u,loc,x,dn3/2                          ! ȥӹ߶
lsel,u,loc,z,0,hc                           ! ȥͲ߶
LESIZE,all, , ,nme, , , , ,1                ! 趨Բͷʷ

vsel,s,loc,x,0,dn3/2                        ! ѡĽӹܼͷ
csys,0                                      ! ֱϵ
vsel,u,loc,x,0,Di                           ! ȥxĲ
vsweep,all                                  ! ɨʷֽӹܼͷxḺ򲿷

lsel,s,radius,,dn1/2-tn1,dn1/2              ! ѡ뾶Ϊӹ1뾶߶
LESIZE,all, , ,nh, , , , ,1                 ! 趨ӹ1ʷ
lsel,r,radius,,dn1/2-tn1                    ! ѡ뾶Ϊӹ1ڰ뾶߶
asll                                        ! ѡ
vsla                                        ! ѡ
aslv                                        ! ѡ
vsweep,all                                  ! ʷֽӹ1
cm,vcon1,volum                              ! ӹ1Ϊvcon1

lsel,s,radius,,dn2/2-tn2,dn2/2              ! ѡ뾶Ϊӹ2뾶߶
LESIZE,all, , ,nh, , , , ,1                 ! 趨ӹ2ʷ
lsel,r,radius,,dn2/2-tn2                    ! ѡ뾶Ϊӹ2ڰ뾶߶
asll                                        ! ѡ
vsla                                        ! ѡ
aslv                                        ! ѡ
vsweep,all                                  ! ʷֽӹ2
cm,vcon2,volum                              ! ӹ2Ϊvcon2

lsel,s,radius,,dn1/2                        ! ѡ뾶Ϊӹ1뾶߶
asll                                        ! ѡ
vsla                                        ! ѡ
cmsel,u,vcon1                               ! ȥvcon1
vsweep,all                                  ! ʷֽӹ1طͷ

lsel,s,radius,,dn2/2                        ! ѡ뾶Ϊӹ2뾶߶
asll                                        ! ѡ
vsla                                        ! ѡ
cmsel,u,vcon2                               ! ȥvcon2
vsweep,all                                  ! ʷֽӹ2طͷ
alls                                        ! ȫѡ
fini                                        ! ˳ǰ
!********  ********
/solu                                       ! 
csys,0                                      ! ֱϵ
asel,s,loc,y,-hc                            ! ѡͲ
da,all,uy                                   ! Լλ
nsla,,1                                     ! ѡڵ
nsel,r,loc,x,0                              ! ѡxΪ0Ľڵ
d,all,ux                                    ! Լxλ
nsla,,1                                     ! ѡڵ
nsel,r,loc,z,0                              ! ѡzΪ0Ľڵ
d,all,uz                                    ! Լzλ
! ********* ѡڱ沢趨Ϊ *************
alls                                        
lsel,s,radius,,Di/2                         ! ѡ뾶ΪͲڰ뾶߶
lsel,u,loc,y,-hc                            ! ȥͲ˲߶
asll                                        ! ѡ
asel,u,loc,y,0                              ! ȥȷ
cm,acon,area                                ! acon

lsel,s,radius,,dn3/2-tn3                    ! ѡ뾶Ϊӹ3ڰ뾶߶
asll                                        ! ѡ
asel,u,loc,y,ky(9)                          ! ȥȷ
asel,u,loc,y,ky(13)                         ! ȥȷ
cmsel,a,acon                                ! acon
cm,acon,area                                ! ٶacon

lsel,s,radius,,dn1/2-tn1                    ! ѡ뾶Ϊӹ1ڰ뾶߶ 
lsel,r,loc,y,ky(49),ky(50)                  ! ѡӹ1߶
asll                                        ! ѡ
cm,acon0,area                               ! acon0
ksel,s,,,49                                 ! ѡؼ
lslk                                        ! ؼѡ
lsel,r,length,,tn1                          ! ѡΪӹ1ȵ߶
asll                                        ! ѡ
cm,acon1,area                               ! acon1
cmsel,s,acon0                               ! ѡacon0
cmsel,u,acon1                               ! acon0ȥacon1
cmsel,a,acon                                ! acon
cm,acon,area                                ! ٶacon

lsel,s,radius,,dn2/2-tn2                    ! ѡ뾶Ϊӹ2ڰ뾶߶ 
lsel,r,loc,y,ky(48),ky(45)                  ! ѡӹ2߶
asll                                        ! ѡ
cm,acon0,area                               ! acon0
ksel,s,,,46                                 ! ѡؼ
lslk                                        ! ؼѡ
lsel,r,length,,tn2                          ! ѡΪӹ1ȵ߶ 
asll                                        ! ѡ
cm,acon1,area                               ! acon1
cmsel,s,acon0                               ! ѡacon0
cmsel,u,acon1                               ! acon0ȥacon1
cmsel,a,acon                                ! acon
cm,acon,area                                ! ٶacon
! ********* ѡӹܶ沢趨Ϊ *************
lsel,s,radius,,dn1/2-tn1                    ! ѡ뾶Ϊӹ1ڰ뾶߶
asll                                        ! ѡ
asel,r,loc,y,ky(59),ky(58)                  ! ѡӹ1
cm,apax1,area                               ! ӹ1apax1

lsel,s,radius,,dn2/2-tn2                    ! ѡ뾶Ϊӹ2ڰ뾶߶
asll                                        ! ѡ
asel,r,loc,y,ky(74),ky(75)                  ! ѡӹ2
cm,apax2,area                               ! ӹ2apax2   

lsel,s,radius,,dn3/2-tn3                    ! ѡ뾶Ϊӹ2ڰ뾶߶
asll                                        ! ѡ
asel,r,loc,y,ky(9)                          ! ѡӹ3
cm,apax3,area                               ! ӹ3apax3
! ********* ѡ沢趨Ϊ *************
lsel,s,radius,,Di/2+tc                      ! ѡ뾶ΪͲ뾶߶                  
lsel,u,loc,y,-hc                            ! ȥͲ˲߶
asll                                        ! ѡ
asel,u,loc,y,0                              ! ȥȷ
cm,aouter,area                              ! aouter

lsel,s,radius,,dn3/2                        ! ѡ뾶Ϊӹ3뾶߶  
lsel,r,loc,y,ky(9)                          ! ѡӹ3˲߶
asll                                        ! ѡ
asel,u,loc,y,ky(9)                          ! ȥȷ
cmsel,a,aouter                              ! aouter
cm,aouter,area                              ! ٶaouter

lsel,s,radius,,dn1/2                        ! ѡ뾶Ϊӹ1뾶߶
lsel,r,loc,y,ky(55),ky(53)                  ! ѡӹ1˲߶
asll                                        ! ѡ
asel,u,loc,y,ky(55),ky(53)                  ! ȥȷ  
cmsel,a,aouter                              ! aouter
cm,aouter,area                              ! ٶaouter  

lsel,s,radius,,dn2/2                        ! ѡ뾶Ϊӹ2뾶߶
lsel,r,loc,y,ky(69),ky(71)                  ! ѡӹ2˲߶
asll                                        ! ѡ
asel,u,loc,y,ky(69),ky(71)                  ! ȥȷ
cmsel,a,aouter                              ! aouter
cm,aouter,area                              ! ٶaouter
! ********* ¶غɲز *************
time,1                                      ! һغɲӦѹ
cmsel,s,apax1                               ! ѡapax1 
sfa,all,1,pres,pax21                        ! Խӹ1ʩӶ˲ƽغ 
cmsel,s,apax2                               ! ѡapax2 
sfa,all,1,pres,pax22                        ! Խӹ2ʩӶ˲ƽغ
cmsel,s,apax3                               ! ѡapax3  
sfa,all,1,pres,pax23                        ! Խӹ3ʩӶ˲ƽغ
cmsel,s,acon                                ! ѡacon
sfa,all,1,pres,p2                           ! ʩѹ
alls                                        ! ȫѡ
solve                                       ! 
save                                        ! 

time,2                                      ! һغɲӦ͹ѹ
alls                                        ! ȫѡ
SFADELE,all,1,pres                          ! ɾغ
cmsel,s,apax1                               ! ѡapax1 
sfa,all,1,pres,pax11                        ! Խӹ1ʩӶ˲ƽغ
cmsel,s,apax2                               ! ѡapax2 
sfa,all,1,pres,pax12                        ! Խӹ2ʩӶ˲ƽغ
cmsel,s,apax3                               ! ѡapax3
sfa,all,1,pres,pax13                        ! Խӹ3ʩӶ˲ƽغ
cmsel,s,aouter                              ! ѡaouter
sfa,all,1,pres,p1                           ! ʩѹ
alls                                        ! ȫѡ 
solve                                       ! 
save                                        ! 
fini                                        ! ˳
! ******    ***********
/post1                                      ! 
LCDEF,1,1                                   ! һغɲΪload case 1
LCDEF,2,2                                   ! ڶغɲΪload case 2 
LCASE,1                                     ! load case 1   
LCOPER,sub,2                                ! ȥload case 2
LCWRITE,3                                   ! Ϊload case 3
PLNSOL,S,INT,0,1                            ! ʾӦͼȡӦǿȷΧ
PATH,a1,2                                   ! 趨·
PPATH,1,48630                               ! 趨·һڵ
PPATH,2,47980                               ! 趨·ڶڵ
PRSECT, ,0                                  ! ȡPLPbQΧжƣ
fini
/post1
! ***** ƣͷ
FTSIZE, 1, 1, 2                             ! 趨ƣλ¼غ
!*  
FP,1,1e6,2e6,5e6,1e7,2e7,5e7                ! ƣSN
FP,7,1e8,1e9,1e10, , ,   
FP,13, , , , , ,
FP,19, ,
FP,21,194,157,127,113,105,99
FP,27,97,96,94, , ,
FP,33, , , , , ,
FP,39, ,
!*  
FL,1,48630,1.0,1.0,1.0,try                  ! ƣͷ
SET,1,last                                  ! һغ
FSNODE,48630,1,1                            ! 㲢洢ƣͷڵĸӦ
SET,2,last                                  ! ڶغ
FSNODE,48630,1,2                            ! 㲢洢ƣͷڵĸӦ
!FE,1,-1                                     ! ǰƣͲ     
FE,1,3.88e5                                 ! 趨¼ѭ
FTCALC,1,48630                              ! ƣ
fini


